Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN ETjiX style file v2.2) 



Limits on the 2.2 contrast ratio of the close orbiting planet 
HD 189733b 

J.R. Barnes 1 * Travis S. Barman 2 , L. Prato 2 , D. Segransan 3 , H.R.A. Jones 1 , C.J. Leigh 4 , 
§ ; A. Collier Cameron 5 and D.J. Pinfleld 1 

i 1 Centre for Astrophysics Research, University of Hertfordshire, Hertfordshire ALIO 9AB. UK. 
• 2 Lowell Observatory, Planetary Research Centre, 1400 West Mars Hill Road, Flagstaff, AZ 86001. USA 
3 Observatoire de Geneve, 51 Chemin des Maillettes, 1290 Sauverny. CH. 
^ , 4 Astrophysics Research Institute, Liverpool John Moores University, Birkenhead CH41 1LD. UK. 
■ 5 SUPA, School of Physics and Astronomy, University of St Andrews, Fife KY16 9SS. UK. 



Accepted 2007 August 27. Received 2007 August 24; in original form 2007 July 20 



ABSTRACT 
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■ We obtained 238 spectra of the close orbiting extrasolar giant planet HD 189733b with 

| resolution R ~ 15, 000 during one night of observations with the near infrared spectrograph, 

NIRSPEC, at the Keck II Telescope. We have searched for planetary absorption signatures in 
the 2.0 |J.m - 2.4 \im. region where H2O and CO are expected to be the dominant atmospheric 
opacities. We employ a phase dependent orbital model and tomographic techniques to search 
for the planetary absorption signatures in the combined stellar and planetary spectra. Because 
potential absorption signatures are hidden in the noise of each single exposure, we use a model 
, list of lines to apply a spectral deconvolution. The resulting mean profile possesses a S/N ratio 

■^J- ■ that is 20 times greater than that found in individual lines. Our spectral timeseries thus yields 

spectral signatures with a mean S/N = 2720. We are unable to detect a planetary signature 
at a contrast ratio of logio(-F p /-F*) = —3.40, with 63.8 per cent confidence. Our findings 
are not consistent with model predictions which nevertheless give a good fit to mid-infrared 
\ observations of HD 189733b. The 1— a result is a factor of 1.7 times less than the predicted 

2.185 |xm planet/star flux ratio of logioX-Fp/i 71 *) ~ —3.16. 
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1 INTRODUCTION 

A factor governing the atmospheric physics of the gas giants in our 
own solar system is the relatively low temperatures of the surround- 
ing environment. While internal heat sources and rotation play a 
role in atmospheric dynamics, much of the structure seen in the 
atmospheres is governed by chemical species which appear in rel- 
atively low abundance (~ 1 per cent by mass) compared with a 
dominant makeup of Hydrogen and Helium. Extra-solar planetary 
atmospheres amenable to study must however differ considerably 
from the solar system paradigm on account of their close proximity 
to their parent star. 

Modern instrumentation is now capable of making ob- 
servations of Close-orbiting Extrasolar Giant Planets (CEGPs) 
in order to learn about the physi cal make up of their at- 
mospheres. Reflected light studies 1 Collier Cameron et al. 19991: 
Charb onneau et al.ll 1999l ; ICollier Cameron et al.ll2002l ; lLeigh et all 
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l2003al lbl) and more recent MOST photometry iRowe et al.l2006h re- 
sults have placed albedo upper limits on the atmospheres of CEGPs 
in the optical. While reflected light studies only search for an atten- 
uated copy of the stellar light, they have the potential to yield im- 
portant information about the chemical makeup of the planetary at- 
mospheres. These upper limits to the planet/star contrast ratio have 
only been able to constrain CEGP atmospheres as having low re- 
flectivity when compared with the solar system gas giants. 

The contrast between the planet and star in the infrared of- 
fers a more favourable regime in which to study CEGPs. An esti- 
mate of the temperature of a CEGP was first reporte d from 24 \im 
second ary eclipse measurements of HD 209458b bv lDeming et al.l 
(2005). While a number of subsequent mid-infrared observations 
of oth er systems have now bee n made with the Spitzer Space Tele- 
scope, Har rington et alj ( 120061) were the first to report measurement 
of the phase dependent flux of a CEGP, implying a world with a sig- 
nificant day-night temperature gradient. Such observations are now 
fuelling new research into the meteorology of extra solar planets 
dCooper & Showmaij|2005ll2006l : lFortnev et alfeOOfj) . To obtain a 
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fuller understanding of the atmospheres of CEGPs it is important to 
study as wide a region of wavelength space as possible on CEGPs 
with a range of parameters. 

The detection of narrow absorption features due to molecular 
species is difficult in the optical because the atmospheres appear so 
dark whereas no high resolution instrumentation in the thermal in- 
frared is widely available. A number of contemp orary instruments 
such as NIRSPEC at K eck dMcLean et alll 19981). IRCS at Sub aru 
dKobavashi et alj200d) and CRIRES at VLT dKaiifl et al.l2004l) are 
capable of making searches for absorption features in the near in- 
frared possible. 

1.1 HD 189733 

HD 189733b was the fir st planet to be disc overed as a transiting 
planet spectroscopically teouchv et ai] |2005). Observations taken 
around the time of transit indicated radial velocity anomalies due to 
the Rossiter-McLaughlin effect, where the transiting planet occults 
part of the stellar disc leading to a deviation from the measured 
Keplerian solution. Further spectrosco pic observations were used 
to determine the spin-orbit alignment ( Winn et al. 2006), indicating 
that the sky projections of th e stellar spin axis an d orbit normal are 
aligned to with 1.4° ± 1.1°. iBakos et al.l {2006) used astrometric 
and radial velocity m easurements to show that HD 189733 is the 
(K1V - K2V) primary (Bou chv et alj|2005f) in a double star system 
with the secondary mid-M dwarf companion orbiting at 216 AU 
from HD 189733 in a plane orthogonal to the orbit of HD 189733b. 

The system pa rameters derived in the discovery paper by 
iBouchv et all [2005) have since been refined through inclusion o f 
multiband photo metric transit obs ervations bv lBakos et al.l (2006). 
More recently, IWinn et alj d2007l) have also used photometric ob- 
servations out of eclipse to determine a stellar rotation period of 
13.4 d. A total of eight transits were used to determine the stellar 
and planetary radii and the photometric mid-transit ephemeris of 
H189733b which are R t = 0.753 ± 0.025 R @ , R p = 1.156 ± 
0.046 Rj, and T c = 2453988.80336(23) + 2.2185733(19) d 
respectively (where the figures in brackets for the ephemeris de- 
note the uncertaint y in the last two d ecimal places). The values 
of 95.79 ± 0.24° dBakos et all l2006h and 95.76 ± 0.29° deter- 
mined for the inclination of the orbit of HD 189733b are in good 
agreement and remove the ambiguity in the planet's mass, yielding 
Mp = 1.13 ± 0.03 Mj. 

Because HD 189733b is one of the closest orbiting planets 
(a = 0.0312 ± 0.0004AU) with a relatively late spectral type par- 
ent star , on e would expec t a rela tively high planet/star contrast ra- 
tio, Fp/F jDeminget al. (2006) detected the strong infrared ther- 
mal emission in the 16 p_m band using the Spitzer space telescope. 
Photometry enabled a clear detection of the eclipse with a depth of 
0.551 ± 0.030 per cent, or F p /F* = 0.0055 ± 0.0003, indicating 
the star to be only 181 times brighter than the planet at 16 [im near 
to e clipse. 

lTinettieTaT] r2007a) modelled transmission spectra, which are 
dominated by H2O and CO opacities, and presented simulations to 
show that these features should be deep enough to produce sig- 
natures which are detectable in t he thermal infrared by present 
and future space based missions. iTinetti et alj {2007b ) have sub- 
sequently used planetary transit depth observations made with the 
Spizter S pace Telescope Infrared Arra y Camera (IRAC) at 3.6 [im, 
5.8 u.m dBeaulieull2007l) and 8 p.m dKnutson et ail 120071) in con- 
juncti on with transmiss ion spectra modelled using H2O opacity 
lists dBarber et ai]|2006t) t o infer the presence o f water in the at- 
mosphere of HD 1 89733b. iGrillmair et aljd2007l) however have re- 



ported measurement of the 7.5 - 14.7 p.m spectrum of HD 189733b 
using the Spizter Space Telescope Infrared Spectrograph (IRS ) 
which does not show the expected decrease dBurrows et alj 2006) 
in flux at the short wavelength end of this range due to incre asing 
H2O opacity. It has been suggested dFortnev & M arlev 2007]) that 
the discrepancy between the IRAC 8 [im flux and IRS flux may re- 
sult from unreliable reduction of the IRS spectra and that the IRAC 
photometry is thus in good agreement with one dimensional atmo- 
spheric models of HD 189733b in the 8 |J.m- 16 |i.m range spanned 
by I RAC flux measureme nts. 

Knuts on et al. U2007l) observed 0.12 per cent flux variation be- 
tween transit and secondary eclipse, with a maximum 16° before 
secondary eclipse. Their 8 p_m Spitzer observations were used to 
derive a simple surface map showing a hot spot offset by 30° East 
of the substellar point. 

In this paper we aim to detect the absorption signature of 
HD 189733b at ~ 2.2 ^.m. The method outlined in Sfjhas the 
potential to determine the plane t/star contrast ratio , F p /F*, if an 
absorption signature is detected (Ba mes et alj2007T) . The detection 
of such opacities forms an important test for the presence of atmo- 
spheric molecular species. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were made with NIRSPEC dMcLean et ail 1 19981) 
at the Keck II Telescope on UT 2006 July 22. A high ca- 
dence timeseries of HD 189733 spectra were recorded using a 
1024x1024 InSb Aladdin-3 array. The rapidly rotating bright A0V 
star, HD 192538, was observed at the start of the night to enable 
identification of telluric features in the spectra. With the NIRSPEC- 
7 blocking filter, a wavelength span of 2.031 1 (im- 2.3809 p_m was 
achieved at a resolution of R ~ 15000. The observations are sum- 
marised in TableQ] 

2.1 Data extraction 

Each raw frame was corrected by subtracting a scaled dark current 
frame. Pixel to pixel variations were removed using flat-field ex- 
posures taken with an internal tungsten reference lamp. In order to 
create a reliable balance frame to remove the pixel sensitivity vari- 
ations, we divided a Gaussian blurred (using a FWHM of 7 pixels) 
version of the master flat field image by the original master flat field 
image. Two flat fields each comprising of 100 median-combined 
frames taken at both the start and end of the night were used since 
a software crash at the end of the first sequence of frames listed in 
TableQ] resulted in an instrument reset. After reboot a small shift 
of 2 pixels in wavelength position of the grating occurred. Exam- 
ination of the balance frames shows that while blemishes and hot 
pixels remain unchanged, a striped or interference pattern present 
in the balance frame taken at the start of the night has shifted in 
the balance frame taken at the end of the night. The interference 
pattern varies from order to order and repeats with a length scale 
of between 10 and 20 pixels. If this pattern is purely dependent on 
the light path, it should largely be removed as a result of flatfield- 
ing. By dividing one balance frame by the other, and measuring the 
residual in each order, we estimate the Hatfield stability limited by 
shifts in the interference pattern after instrument setup changes (2 
pixel shift) through the night to be 1.0 per cent. The Hatfield at the 
start of observations was used for the first block of observations 
which only comprised 8 sets of co-added frames (TableQ] see cap- 
tion), while the Hatfield taken at the end of observations was used in 
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Table 1. NIRSPEC/Keck II observations of HD 18733 and HD 192538 for UT 2006 July 22. Cloud cover was highly variable throughout the night. Groups of 
15 spectra from the first HD 189733 block were coadded post-observations to give the same effective exposure time (60 sees) as in the second block, resulting 
in an effective total of 238 frames. As a result of variable cloud cover, we were unable to extract 1 3 of the HD 1 89733 spectra, leaving 225 usable frames. 



Object 


UT start of 
first frame 


UT start of 
last frame 


Time per Number of 
exposure [sees] coadds per frame 


Number of 
observations 


Comments 


UT 2006 July 22 


HD 192538 
HD 189733 
HD 189733 


05:48:46 
06:10:58 
07:04:27 


05:59:18 
06:43:52 
14:08:0922 


60 1 
4 1 
4 15 


8 

125 

230 


A0V template 

Software crash at end of sequence 



conjunction with the remainder of observations after the instrument 
was reset. 

The worst cosmic ray events were removed at the pre- 
extraction stage using the Starlink FIGARO routine BCLEAN 
dShortridgel 19 931 The recorded frames comprised 4 sec exposures 
at the start of the night before the software crash (see TableQJ and 
15 x 4 sees = 60 sees exposures thereafter. We therefore coadded 
groups of 15 exposures from the first set of frames to ensure that 
all spectra were exposed for the same 60 sees. 

A number of sky lines were present in the four shortest 
wavelength orders while the two reddest orders contained only one 
very weak line. In order to maximise the stability of observations 
made with the instrument, we decided against using a nodding 
ABBA sequence. We were thus unable to use subsequent A-B 
and B-A differences to reject sky lines during t he extract i on. W e 
followed a similar method to that detailed in lBarnes et alj d2007l) , 
finding that the Starlink echelle data reduction package, ECHOMOP 
(Mills 119940 . rejects all but one strong sky line, at ~ 21501.5 A 
through iterative fitting of a third degree polynomial. Bad pixels 
on the detector close to the location of this strong sky line meant 
that a small region encompassing the line was excluded from 
further analysis. The spectra were extracted using ECHOMOP's 
im plemen t ation of the optimal extraction algorithm developed 
by iHornej qJ986). ECHOMOP propagates error information based 
on photon statistics and readout noise throughout the extraction 
process. 



2.2 Wavelength calibration 

Calibration frames using a combination of Ar, Ne, Xe and Kr 
lamps were taken during the observations although the line iden- 
tifications are p oor in some orders . We therefore adopted the ap- 
proach used in iBarnes et al. I d2007t) where telluric features in the 
observed spectrum of the standard star were used. A spectrum 
generated from a HITRAN line list teothman et ail 2005) enabled 
identification of the corresponding features with known/measurable 
wavelengths. Between 10 and 19 lines were used in each or- 
der and a cubic polynomial was used to give fits with a rms 
scatter equal to 1.44 x 10~ 6 tim (0.0144 A). The mean wave- 
length increment of the six calibrated orders (NIRSPEC orders 
32-37) is 3.23 x 10~ 5 [im (0.323 A). At the 2.185 \im cen- 
troidal wavelength position of our deconvolved lines (see i]3.3t 
the rms scatter in the fit corresponds to 0.045 pixel elements or 
0.00099 resolution elements. The wavelength ranges for each of the 
six orders are: 2.031 -2.061 [im (order 37), 2.087-2.118 (im (or- 



der 36), 2.146-2.178 ^m (order 35), 2.209-2.242 tim (order 34), 
2.275 - 2.309 \im (order 33), 2.346 - 2.381 urn (order 32). 



3 SEARCHING FOR A PLANETARY SIGNATURE 

3.1 Removal of stellar and telluric spectra 

Before searching for a planetary signal, it is important to remove 
the stellar and telluric lines from our spectra. This is achieved by 
subtracting a shifted and scaled template spectrum that has been 
created by coadding all the frames from our night of observations. 
This procedure is valid because the planetary signal is not station- 
ary during the observations. The phases of the observations made 
here mean that any planetary absorption signature will be Doppler 
shifted from +60kms _1 to -80kms _1 during the planet's or- 
bit about HD 189733. The template spectrum was fitted to each 
HD 189733 spectrum in turn by taking derivatives of the spectra 
and using splines to calculate the scale factor at points across the 
spectra. This process can account for those lines that behave inde- 
pende ntly over the night and is described in Coll ier Cameron et al.l 
hoo2h (Appendix A). For example, all telluric lines do not neces- 
sarily vary in strength by the same factor at any given time while 
the strength of all telluric lines may vary relative to the stellar lines. 
Fig. 1 (tick marks at bottom of each panel) shows the positions of 
stellar lines (f or further details s ee £|3 ,31 > from the Vienna Atomic 
Line database iKup kaet alJl999h .A comparison with the observed 
template spectrum above the tick marks shows that the spectra are 
dominated by the telluric lines. 

While this template subtraction procedure is effective at re- 
moving the stellar and telluric lines and preserving a possible plan- 
etary signal, some residuals remain in the spectra. These residuals 
often appear in blocks of spectra and become evident when each 
echelle order is displayed in a greyscale image. The large variabil- 
ity in S/N ratio of the current spectra ( S]3.3I >. due to highly variable 
cloud cover, prompted us to investigate removal of the blocks of 
spectra containing the fewest counts. In Sj4]we report our findings 
for the case where 42 per cent of the spectra were rejected prior to 
further analysis. 

Remaining systematics w ere removed using principal c ompo- 
nents analysis as described bv lCollier Cameron et al. 1 20021) (Ap- 
pendix B); and discussed at length in IBarnes et alj < |2007|) . A bal- 
ance must be struck between removal of systematics and removal 
of photons which would reduce our detection threshold. Only resid- 
uals at a fixed wavelength position are removed if the first few prin- 
cipal components are subtracted from the data, ensuring a moving 
planetary signal is not significantly attenuated. We subdivided each 



Co 



3 




Model Planet 




Observed spectrum 



■iiBiiiiiiiiiiiiiiiiiiii i iiiiiiiii iiiiiiiiiii nun i mi nun ii in i ii iiii m i mill ii inn i mill urn 11 n u i mi in mi i inn i u i 1 1 urn i m i n inn n n i 



H„0 



1 1 1 iii an iiiiiiiiiii ii iiii 1 1 iiiiiiiiiii mil ii ii ii i mi 



2.04x 1 



2.06x1 CT 



2.08x10 4 



2.1X1CT 



2.12x10 4 



2.14x1 0* 



2.1 6x1 4 



Wavelength [A] 



2.1 8x1 CT 



2.2x10 4 



Model Planet 




CO 



HF 



i !■ ii 1 1 1 ii ■! ii 1 1 ii i in ii mi 



mini i ii i i in ii ii 

iiiiiiiii i i i i i i inn i i 



Observed spectrum 




H 2 



2.22x10* 



2.24x10 4 



2.26x10* 



2.28x10* 2.3x10* 
Wavelength [A] 



in i 
2.32x10 4 




2.34x1 0* 



2.36x10* 



2.38x10* 



Figure 1. Model planet (top in each panel) spectrum and observed spectrum (bottom in each panel) for NIRSPEC orders 37, 36, 35, 34, 33 & 32 (top to bottom, left to right). Each 
spectrum spans 2.03 ^m - 2.38 ^m at a resolution of R ~ 15000. The model planetary spectrum has been degraded to the same instrumental resolution and broadened with a rotation 
profile of width 2.53 km s _1 in accordance with a tidally locked planet with the empirically determined parameters of HD 189733b. The gaps between spectra show the regions not 
covered by the instrumental setup (see §2.2). The planetary spectrum is shown in the absence of any stellar host and in reality would be weaker by a factor equal to the planet/star 
contrast ratio, F p /F*. The tick marks at the bottom of each panel represent the theoretical positions of Fe and Si stellar lines as given by VALD for a T — 5000 K spectrum. The tick 
marks in the middle of each panel mark the positions of the strongest expected planetary opacities from our model (labelled accordingly from bottom to top as H2O, CO and HF) 
demonstrating H2O to be the dominant species. The HF and CO opacities begin at 2.2895 fim and 2.2935 fim respectively. 
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of the six orders into three subsections and subtracted the strongest 
components (between 2 and 5) for each section. 



3.2 Atmospheric Model 

We attempt to make optimal use of the information content of our 
spectra by considering all the planetary absorption lines contained 
within the ~ 1645 A wavelength span of our data. We therefore re- 
quire a mod el that describ e s the p lanetary atmosphere as accurately 
as possible. IBarnes et alj 12007) discuss the nature of the irradi- 
ated model atmosphere that includes day-night time temperature 
gradients on the planet. For a ful l descr i ption of the model opaci- 
ties and setup see Ferguson et alj ( l2005l)jBarman etHEI bOOll) and 



iBarman et at] d2005h BHA05. For HD 189733b, we model an at- 
mosphere with solar metallicity for a temperature, T = 1250 K and 
surface gravity, log g = 1.33 ms~ . 

In [Barnes et all J2007h . the modelled planet, HD 75289b, 
was not as close to the GOV parent star (a — 0.0482 AU) as 
HD 189733b is to it's K2V host (a = 0.0312 AU). The opacities 
that remain after rainout of condensed species has occurred (such as 
TiO and VO) appear stronger for HD 1 89733b when compared with 



liuand VU) appear stronger tor tiU 189/ Jib when compared with 
HD 75289b ( IBarnes et al. 120071) and are due to two main species, 



namely H2O and CO, in the wavelength region of our data. Fig. 12.21 
shows the wavelength range of our data where the tick marks rep- 
resent strongest H2O, CO and HF opacities. The detection method 
used to search for a planetary signature ( S]3.3l and S]3.4t is model 
dependent and may be affected by a number of factors. Missing or 
incorrect line opacities in the model are potential problems; how- 
ever, the line opacities used in the model are known to be fairly reli- 
able for similar tempe ratu res and pressures f ound in brown dwarfs 
dFreedman et alj2007l) . In lBarnes et alj d2007l) . we investigated the 
effects of a mismatch in line opacity strengths and line positions 
for a similar instrumental setup to that used for the observations 
in this paper. We found that the simulations still enabled us to de- 
tect a planetary signal with 99.9 per cent confidence (a 1.58 times 
drop in planet/star ratio) for the case where the model line list line 
strengths used in deconvolution were modified by 50 per cent or for 
the case where the 15 percent of the line opacity wavelengths were 
randomised. The H 2O line lists used in our model are taken from 
iBarber et alj d2006l) and represent the most complete and accurate 
water line lists to date. 



3.3 Deconvolution 

A line list (comprising opacity wavelength positions and line 
depths) based on the model atmosphere described in i]3.2l is used as 
a deconvolution template for planetary absorption features which 
are hidden within each HD 189733b spectrum. The least squares 
deconvolution procedure yields a single mean absorption line for 
each HD 189733b spectrum which when scaled with the line list 
gives the best fit to the observed spectrum. In this way the effective 
S/N ratio of the mean deconvolved absorption line is much higher 
than that of individual lines, potentially enabling a weak signal to 
be extracted from the data. The method of least squares deconvo- 
lution was first describ ed by iDonati et alj dl997t) . Our implemen- 
tation of the algorithm ( IBarnes et al.ll 19981) propagates errors from 
the input spe ctra and has been used in reflected light searches in 
the opti cal bv lCollier Cameron et ail dl999ll2002h and lLeigh et~ai1 
(2003a b). A detailed description o f the use of this algo rithm with 
near infrared spectra is described in Barne s et al. (2007). 

Of the 225 usable spectra, the mean S/N ratio, measured from 



residual spectra after removal of the template, was 136 ± 85. The 
minimum and maximum S/N ratios were 27 and 382 respectively 
with only 10 per cent of the spectra possessing S/N > 250. We 
assess the consequences of this highly variable S/N in section fj4] 
Since the spectra are dominated by very strong telluric lines in the 
2.0445 - 2.0610 11m region of order 37 and the 2.3650 - 2.3807 nm 
region of order 32 (Fig. I2.2I >. we did not use these regions in the 
deconvolution procedure. The gain in S/N of 20 for the mean de- 
convolved spectrum yielded S/N ratios of 2751 ± 1772. The full 
range of S/N ratios in the deconvolved spectral profiles is however 
540-7640, again, reflecting the highly variable cloud conditions 
throughout the night of observations. 



3.4 Matched filter 

We model the time-dependent changes in Doppler shift and 
brightness using a matched filter analysis described in detail in 
Appendix D of lCollier Cameron et alj J2002I). Mo difications of this 
implementation are described in lBarnes et ail d2007l) . but we reiter- 
ate the basics of the procedure here for clarity. 

We model a travelling Gaussian which traces out the motion of 
HD 189733b in an assumed circular orbit about its host star. The ra- 
dial velocity motion of the planet as a function of phase thus traces 
out a sinusoid. We assume the planet to be tidally locked and use a 
model describing the variation in brightness of an irradiated planet 
that possesses a hot inner face with littl e redistribution of heat. We 
assume, as found in lBarnes et alj d2007l) , that this vari ation closely 
resembles a Venus-like phase function dffiltorJll992l) . Discussion 
of the true nature of the phase function is discussed in more detail 
in [15] in light of recent observations, although we note here that a 
Venus-like phase function agrees with these findings to within 5 
per cent over the range of phases of our observations. We asses the 
relative probabilities of the \ 2 fits to the data by varying the peak 
contrast ratio, eo(A), and the velocity amplitude of the orbital mo- 
tion, K p . The peak contrast ratio occurs at secondary eclipse (i.e. 
phase <f> = 0.5) in the case of a Venus-like phase function. In the 
phase range (0.44 < (f> < 0.59) of our data, we are sampling the 
region of the sinusoid where, to first order, the radial velocity mo- 
tion approximates a straight line. The gradient of this line, which 
is fitted by a sinusoid in our model, thus enables an estimate of 
the velocity amplitude, K p , to be made. The improvement in \ 2 ls 
estimated for combinations of these parameters and is normalised 
relative to the best-fitting model. To calibrate any candidate detec- 
tion, we construct a simulated planet signal of known eo(A), with 
A = 2.185 \im that is added to the extracted spectra prior to analy- 
sis. By ensuring the fake planet is recovered at the correct contrast 
ratio (hereafter, simply eo or logio eo) by our procedures, we can be 
confident of a detection in the presence of a genuine planet signal. 

The significance of the result is assessed using bootstrap sta- 
tistical procedures based on the random re-ordering of the data in a 
way that scrambles phase s while preserving the effects of any cor- 
related systematic errors dCollier Cameron et alj [2002"). The order 
of the observations is randomised in a set of 3000 trials which will 
scramble any true planetary signal while enabling the data to retain 
the ability to produce spurious detections through the chance align- 
ment of systematic errors. The least squares estimate of logio eo 
and associated \ 2 as a function of K p enable us to plot 68.4, 95.4, 
99.0 and 99.9 per cent bootstrap limits on the strength of the plan- 
etary signal. 
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Figure 2. Left: Phased timeseries of the deconvolved residual spectra of HD 189733. The dashed sinusoidal curve represents the expected motion of a 
planetary signal based on the empirically determined velocity amplitude, K v = 152.8 ± 2.0 km s _1 . Since the S/N ratio is variable, for plotting purposes, all 
deconvolved profiles have been normalised to the same noise level to optimise the visible information content of the plot. Right: Relative probability \ 2 map 
of planet-star flux ratio logio £0 vs K p . The first 2 to 5 principal components have been subtracted from the timeseries spectra. The greyscale in the probability 
maps represents the probability relative to the best fitting model in the range 100 km s" 1 < K p ^ 200 km s" 1 in the range white for to black for 1 . A 
candidate feature with K p = 52 km s _1 resulting from systematics (see main text) has resulted in the large dark region in the left of the plot. Plotted are the 
68.3, 95.4, 99 and 99.9 per cent (bottom to top) confidence levels for a planetary signal detection. The dashed vertical line represents the velocity amplitude, 
K p = 152.8 km s _1 . No planetary signal is detected at this velocity. Top: Analysis using all 225 spectra. Bottom: Analysis after removal of blocks of the 
lowest S/N ratio spectra. Using the remaining 131 deconvolved spectra results in the same confidence levels appearing at lower contrast ratios. 
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4 RESULTS 

In Fig. [2] (left), we plot the deconvolved spectral timeseries. The 
dashed line is plotted as a guide for the reader and represents the ve- 
locity path of a planet with the same parameters as determined em- 
pirically for HD 189733b. We can obtain an estimate of the velocity 
amplitude of the planet by taking a = 0.0312 ± 0.0004 AU, P = 
2.2185733 ± 0.0000019 d and i = 85.76 ± 0.29°, from which we 
find v p — 152.58 ± 1.96 km s" 1 . We discuss additional factors 
that could augment this uncertainty in Sj5] Fig. [2] (right) presents 
the result of our matched Gaussian filter analysis. Again the dashed 
line represents the velocity at which HD 189733b should appear. 
The darkest regions in the map represent the greatest improve- 
ment in x 2 when fitting the model described in 93.31 In the op - 
tical studie s carrie d out by ICollier Cameron et alj dl999l 120021) : 
Leig h et al. I ( 12003 M . the reflected light signal was simply an at- 
tenuated copy of the stellar spectrum. The central portion of the de- 
convolved timeseries was excluded during the matched filter mod- 
elling to ensure that no spurious signal due to remaining residuals 
of an inadequately subtracted stellar profile were included. Here 
however, the planetary spectrum is not expected to be a copy of 
the stellar spectrum and any line opacities should be randomly 
distributed with respect to stellar opacities in the spectrum. How- 
ever, since HD 189733b is a transiting planet, when applying our 
matched filter analysis, we exclude t he po rtion of spectrum when 
the planet is not visible. IWinn et alj fe007h estimate a transit time 
of 1.827 ± 0.029 hrs indicating that the planet is not visible for 
phases 0.483 < <f> < 0.517. 

The feature at K p ~ 52 km s _1 is the result of remaining sys- 
tematics in the deconvolved timeseries. To enhance the greyscale in 
the region of the x 2 ma P where we know HD 189733b should ap- 
pear, as outlined above, we excluded K p < 100 km s _1 when com- 
puting the greyscale, resulting in a dark block in the map around the 
K p ~ 52 km s _1 feature. We do not detect the planet with 1-cr and 
2-cr confidence at levels of logio eo = log w (F p /F, ) = —3.26 
and —2.72 respectively. The equivalent star/planet contrasts are 
F,/F p = 1820 and 525 respectively. 

Fig.|2](top left) shows clear systematic features, which are not 
subtracted, even after removing the first 2-5 principal components 
as outlined in 93.11 This is especially noticeable for spectra taken 
in the phase range 0.52 < c6 < 0.54. We attribute these residuals to 
variation in telluric strengths throughout the night and as a result of 
the highly variable cloud cover. Since the observations were made 
at a spectral resolution of R ~ 15, 000, it is not possible to effec- 
tively scale the template spectrum in all regions, especially where 
those telluric features that may vary in strength independently are 
blended. 

To investigate the effect of the highly variable S/N, we re- 
jected frames with the lowest S/N ratios. A threshold level that re- 
moved the main blocks of lower S/N spectra ( S/N < 97) was cho- 
sen, leaving 131 spectra. The deconvolved profiles (Fig. [2] bottom 
left) possess S/N = 3800 ± 1480. The matched filter analysis per- 
formed on these spectra alone enables us to place 1-cr and 2-cr con- 
fidence at levels of logio eo = —3.40 and —2.88 on a non-detection 
of HD 189733b (Fig. [2] bottom right) . The equivalent star/planet 
contrasts are F,/F p = 2512 and 759 respectively. We are nev- 
ertheless able to improve the sensitivity of our result, although at 
the 1-cr level, the results are indistinguishable. Despite using fewer 
spectra, we believe that an increased sensitivity is due to the re- 
duced range of S/N values of our spectra. The systematic effects 
introduced by unresolved telluric features, which cannot easily be 
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Figure 3. Model planet/star flux ratio for the HD 189733 system. A Spitzer 
eclipse depth measurement u sing the Infrared Array Camera (IRAC) is plot- 
ted for 8 \xm (Knutson et al. (2007), K07). Similar measure ments using the 
Spizte r Infrared Spectrogra ph (IRS) are plotted fo r 11.1 |jm dGrillmair et alj 
120071) . G07) and 16 yon iDeming et all <2006l) , DHSP06). A horizontal 
bar with vertical down-pointing arrow indicates our 2. 1 85 \jm 1 — cr upper 
limit. The width of the horizontal bar represents the wavelength range of 
our data. The single point plotted in the inset at logio ^0 = —3.16 repre- 
sents the model mean flux over this range. The upper limit is at a level of 
logio £0 = 0.24 lower than the model (i.e. eo(model)/eo (observed) = 1.7). 



corrected for dBailev et alj|2007h . are thus a critical factor in deter- 
mining the sensitivity of our technique. 



5 DISCUSSION 

The near infrared planet/star contrast ratio, F p /F t , offers a much 
more favourable regime over optical wavelengths in which to carry 
out a search for a planetary signal. The technique, in addition to 
offering a chance to determine the contrast ratio, relies upon our 
ability to detect absorption in the planetary atmosphere. The con- 
trast limit placed by this study is intriguing because we do not de- 
tect the planet at a level of F p /F, that is 1.7 times less than the 
model prediction; and despite the observation that ou r model is in 
agreement with mid-infrared Spitzer observatio ns dPeming et al.1 
120061: iGrillmair et ai]|2007l ; iKnutson et al]|2007l) . Fig. [3] is a plot 
of the theoretical F p /F* from our model indicating both the mid- 
infrared Spitzer measurements and the upper limi t determined 
above. For comparison with the Bar nes et al.1 J2007t) simulations, 
a non-detection at the reported level may arise if 55 per cent of 
line strengths are mis-matched or 16 per cent of the lines possess 
incorrect wavelengths (see 93. 21 . 

Apart from line position and strength mismatches between 
model and observed spectra, another possibility that would result in 
a non-optimal searc h for the planetary s ignature is use of an incor- 
rect phase function. iBarnes et al.1 d2007h found a Venus-like phase 
function for the HD 75289 models over the phases of observation 
consid ered. The overall variation in flux found by IKnutson et al.1 
(2007) indicates that a Venus-like phase function is less appropri- 
ate for HD 189733b, indicating that heat is more effectively trans- 
ported in its atmosphere. The lower contrast between night and day 
yields a temperature range of 973 - 1212 K. Over the small range of 
our observations however, we find that the variation in flux from the 
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maximum as observed bv lKnutson et alj d2007h is consistent with a 
Venus-like phase function to within 5 per cent. 

lHarrington et al] d2006h measured the phase dependent flux of 
v And b, finding variation that was well fitted by the models of 
BHA05, implying poor transfer of t hermal energy. On th e other 
hand, the more recent observations bv lKnutson et al.l |2007) found 
that a maximum in flux for the HD 189733 system occurred at a 
phase angle of 16° before mid-eclipse of the planet. They derived 
a surface map of HD 189733b and found the hottest part of the at- 
mosphere is offset 30° East of the substellar point. We adjusted our 
Venus-like phase function by 16° to ensure a maximum occurred at 
the same phase (i.e. <f) = 0.456) rather than at phase <f> = 0.5. This 
did not result in the appearance of a candidate planetary feature in 
our x 2 ma P however. The resulting increase of the l-cr confidence 
level from log eo = -3.4 to log to = -3.25 is equivalent to a reduc- 
tion in sensitivity by a factor of 1.4. This is probably the result of 
shifting the time at which the planetary flux maximum occurs to 
the phase where observations began. As can be seen in 

Fig.|2](bot- 

tom left), only one observation is made before maximum when the 
shifted phase function is used. 

Nevertheless, the possibility that the observed phase function 
peak may be a functi o n of w avel ength is supported by simulations . 
ICooper & Showman! d2005l) and ICooper & Showman! d2006t) car- 
ried out three-dimensional atmospheric circulation simulations for 
HD 209458b and found that the circulation pattern is a function 
of atmospheric pressure or depth. In the lower pressure regime of 
the outer atmosphere, the radiative timescale is much shorter than 
advective timescales such that absorbed stellar heat is quickly re- 
radiated. Lower in the atmosphere, at higher pressures, advective 
timescales become important and carry a bsorbed energy in an east- 
erly direction. Based on these models. [Fortnev et al.1 d2006l) has 
shown that the location of the phase function peak is indeed ex- 
pected to be a function of wavelength in the case of HD 209458b. 
The time at which the phase function peak occurs before secondary 
eclipse is inver sely proportional to wav elen gth, in agreement wit h 
the findings of lHarrington et al. I d2006l) and Knutson et al. I d2007h . 
This implies that shorter wavelengths probe deeper layers of the 
planetary atmosphere. Wind speeds of order 6.6 — 4.3 km s _1 at 
10 mbar to 1 bar in an easterly direction are high, and while not 
sufficient to significantly shift the apparent velocity amplitude of 
the planet (e.g. see Fig. [2j, they introduce an additional factor into 
the error budget of the estimated 1.96 km s _1 uncertainty ([Q. An 
easterly wind flow would act to reduce the apparent velocity am- 
plitude, but since most of the planetary signal is expected to come 
from the hotspot which would be travelling transversely to the ob- 
server at the observed phases, any induced perturbation is likely to 
be a second order effect. 

Although iFortnev et al.l d2006l) do not consider wavelengths 
short of 3.6 ]xm it seems reasonable to expect that our ~ 2.2 iim 
observations also probe deeper layers of the atmosphere. If the 
temperature-pressure gradient is lower or even flat in this regime 
(i.e. lOOmbar - 1 bar), as indicated for HD 209458b in BHA05, it 
may be that absorption lines are suppressed and the dayside spec- 
trum approaches that of a black body. Also, the reasonable agree- 
ment between the hotter no-redistribution model and the Spitzer 
observations (that emerge from layers above the AT-band pho- 
tosphere), may be an indication of a deeper temperature struc- 
ture cooled by energy redistribution to the nig ht side, similar to 
what h as been predicted for HD 209458b bv lCooper & Showmanl 
(2005). Constraints on the depth-dependent redistribution in the at- 
mosphere of HD 189733b will be explored in more detail by Bar- 
man et al. (in prep.). 



6 CONCLUSION 

The methods presented here provide a powerful technique with 
which to detect a planetary absorption signature and the l-a 
logio eo = -3.40 limit is extremely interesting. We have outlined 
a number of possibilities for the lack of detection of a planetary 
signal. Our data suggest that ideal observing conditions would en- 
able us to obtain spectra with S/N ratios as high as ~ 380. With 
the same NIRSPEC instrumental setup at Keck II, a total of 230 
frames with S/N = 300 taken over the same orbital phase range 
would enable us to push our 1 — a and 2 — a limits down to levels 
of logio eo = -3.90 & 3.40 respectively. Also, the high variability of 
the stellar flux relative to the telluric lines (due to highly variable 
cloud cover) seen in our present data would be reduced in more 
ideal conditions. This would lead to fewer systematics after tem- 
plate spectrum removal enabling us to further lower our confidence 
level estimates for eo. Additionally, observations at higher spectral 
resolution than currently af forded by the presen t data would be de- 
sirable, as demonstrated in Barn es et al. I d2007l) . and would lead to 
better telluric subtraction. 

An independent, near infrared, photometric detection would 
enable us to place stronger limits on the reliability of the line list pa- 
rameters or the nature of the atmosphere itself. For instance, when 
combined with spectral observations, as described above, it would 
be possible to rule out or place limits on the presence of absorption 
features in the corresponding wavelength region. 
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